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l^j | Abstract 

■ Several explicit examples of quasi exactly solvable 'discrete' quantum mechanical 

Hamiltonians are derived by deforming the well-known exactly solvable Hamiltonians of 
one degree of freedom. These are difference analogues of the well-known quasi exactly 
t^- ■ solvable systems, the harmonic oscillator (with/without the centrifugal potential) de- 

formed by a sextic potential and the 1/ sin 2 x potential deformed by a cos 2x potential. 
They have a finite number of exactly calculable eigenvalues and eigenfunctions. 

: 
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1 Introduction 

X' 

Exactly solvable and Quasi Exactly Solvable (QES) quantum mechanical systems have played 
a very important role in modern physics. The former, the exactly solvable systems, are quite 
well-known. In the Schrodinger picture, if all the eigenvalues and corresponding eigenfunc- 
tions are known, the system is exactly solvable. Plenty of such systems are known, for exam- 
ple, the Poschl-Teller and the Morse potential on top of the best-known harmonic oscillator 
and the coulomb potential [1] for degree one cases and the Calogero-Sutherland systems 
[21 [3j H] for many degrees of freedom cases. Recently the exact Heisenberg operator solutions 
and the corresponding annihilation-creation operators are constructed for most of the degree 
one exactly solvable quantum mechanics [5] and for the multi-particle Calogero systems [6j. 
The notion of exact solvability was generalised to the so-called 'discrete' quantum mechan- 
ics, in which the Schrodinger equation is a difference equation in stead of differential. The 



difference analogues of the Calogero- Sutherland systems were constructed by Ruijsenaars- 
Schneider-van Diejen [7J[8]. The difference equation analogues of the equations determining 
the Hermite, Laguerre and Jacobi polynomials were derived within the Askey-scheme of hy- 
pergeometric orthogonal polynomials [9] [10] . Later they were reformulated as Hamiltonian 
dynamics with shape-invariance [IT] by Odake-Sasaki [T2] . 

In contrast, the latter, Quasi Exactly Solvable (QES) systems have a short history and 
less known. If a finite number of exact eigenvalues and eigenf unctions are known, the system 
is QES [13]. Since the number of exactly solvable states can be chosen as large as wanted, a 
QES system could be used as a good alternative to an exactly solvable system. Several one 
degree of freedom QES systems are listed in [13], |Tj] and multi-particle QES systems were 
first constructed by Sasaki- Takasaki [15] as deformation of Inozemtsev- Calogero- Sutherland 
systems, which was followed by [IB] . 

In the present paper we derive several explicit examples of QES difference equations 
as deformation of exactly solvable 'discrete' quantum mechanics [12] . They are differ- 
ence analogues of the well-known quasi exactly solvable systems, the harmonic oscillator 
(with/without the centrifugal potential) deformed by a sextic potential and the l/sin 2 a; 
potential deformed by a cos 2x potential. 

This paper is organised as follows. In the next section, the deformation method to 
obtain a QES from an exactly solvable system is explained in some detail by taking two well- 
known examples of the ordinary quantum mechanical QES systems. Then two corresponding 
difference equation analogues are derived. Section [3] provides three more explicit examples. 
The final section is for a summary and comments. 

2 Quasi Exactly Solvable Deformation 

There are many different ways of deriving QES Hamiltonians for ordinary quantum me- 
chanics [T3l [H] . However, to the best of our knowledge, a very limited number of explicit 
examples of QES difference equations are known in connection with U q (sl(2)) [17] . In these 
examples, quantum wavefunctions are related to those defined on discrete lattice points only. 
In the present paper we present several explicit examples of QES 'discrete' quantum mechan- 
ical Hamiltonians of one degree of freedom, whose wavefunctions are continuous functions of 
x as in the ordinary quantum mechanics. They are obtained by deforming exactly solvable 
'discrete' quantum mechanical Hamiltonians [12], which have the Askey-scheme of hyper ge- 
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ometric orthogonal polynomials [9l [10] as part of the eigenf unctions; the Meixner-Pollaczek, 
continuous Hahn, continuous dual Hahn, Wilson and Askey- Wilson polynomials. This de- 
formation method was first applied by Sasaki and Takasaki [15] to derive multi-particle QES 
based on the Inozemtsev models. 

For illustrative purposes, we will explain the deformation method for the two well-known 
examples of degree one QES systems in ordinary quantum mechanics in the next subsection. 
These examples are the sextic (x 6 ) potential added to the harmonic oscillator (x 2 ) potential, 
and another a cos 2x potential added to the exactly solvable 1 / sin 2 x potential. The same 
method is applied in subsection 12.21 to derive the first two examples of QES Hamiltonian in 
'discrete' quantum mechanics corresponding to the the sextic potential deformation. The 
rest of the examples are given in section [3j 

2.1 Ordinary Quantum Mechanics 

2.1.1 Harmonic Oscillator Deformed by Sextic Potential 

A best-known example of QES Hamiltonian, the harmonic oscillator plus a sextic (x 6 ) po- 
tential is given succinctly by 

W=-- J - J + — +-— + a M {x), a M (x) = -2aMx 2 , MeN. (2.1) 
dx z \ dx J dx z 

Here we tentatively call the last term of the above Hamiltonian, the compensation 

term. Throughout this paper we adopt the unit system 2m = h = 1. The real prepotential 
W is a deformation of that for the harmonic oscillator Wq = —bx 2 /2, with a being the 



deformation parameter: 



def 0- 4 T -., O 4 ^2 



W = W{x) = --x 4 + W = --x 4 - -x\ a, b > 0. (2.2) 

By the similarity transformation in terms of the pseudo ground state wavefunction <fio(x) = 
e w ( x \ we obtain 

~ def , i „ , d 2 dW(x) d , . 9 . . 

H = 0Q 1 oHo0 o = - -2— -2a>ix 2 , (2.3) 

lAltAJ \JUiAj \XitXj 

d 2 d 
= -—r+(2ax 3 + 2bx) — -2aMx 2 . (2.4) 

dx 1 dx 

In the absence of the compensation term a^fi), <fio(x) is actually a ground state wave- 
function, therefore it has no node and it is square integrable. Another characterisation of 
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the pseudo ground state wavefunction 0o is that it is annihilated by the operator A which 
factorises the main part of the Hamiltonian (12.11) : 

AM*) =0, A^-± + ^- t A^±+™2£, (2.5) 

ax ax ax ax 

H = A^A~2aMx 2 . (2.6) 

The action of the Hamiltonian Ti, (12. 3p on monomials of x reads 

( -n(n-l)x n - 2 + 2nbx n + 2a(n-M)x n+2 , n < M - 2, 
H x n = < _ (2 7) 

[ -M(.M - l)x M ~ 2 + 2Mbx M , n = M. 

Since the parity is conserved, it is now obvious that H keeps the polynomial space Vm 
invariant, 

HV M C V M , (2.8) 

v def f Span [l, x 2 , . . . , x 2k , . . . , x M ] , M : even, g ^ 
I Span [x,x 3 , . . . ,x 2fc+1 , . . . ,x M ] , A< : odd, 

and that Ti. is a tri-diagonal matrix. Thus we can obtain a finite number of exact eigenvalues 
and eigenf unctions of the sextic potential Hamiltonian (12.11) in the form: 



n<\> = £<f>, 4>{x) = o (x)Pm(x), p m ev m ,^ np M = £P M , (2.10) 

by the diagonalisation of a finite dimensional Hamiltonian matrix 7i (12.31) with 

{ M/2 + 1, .M:even, 
dimV ^ = ( (M + l)/2, ^:odd. (2 - U) 

Since 7i is obviously hermitian (or self-adjoint), all the eigenvalues are real and eigenfunctions 
belonging to different eigenvalues are orthogonal with each other. In other words, two 
polynomial solutions Pm and P' M are orthogonal with respect to the weight function 4>q(x). 
The square integrability of all the eigenfunctions of the above form (I2.10p (j) 2 (x)dx < oo is 
obvious. The true ground state wave function has the form fl 2 . 1 [) with the lowest eigenvalue, 
say Sq and it has no node due to the oscillation theorem. 

2.1.2 1/ sin 2 x Potential Deformed by cos2x Potential 

Another well-known example of quasi exactly solvable system is the exactly solvable 1 / sin 2 x 
potential (Wo == g log sin x) deformed by a cos2x potential. The Hamiltonian has the same 



form as (12.11) with only the prepotential W(x) and the compensation term cx,m{ x ) different: 

H = ~^+(^X + ^ + <* M (x), = 4aM sin 2 x, M.EN, (2.12) 

dx A \ ax J dx z 

W d = - cos2x + Wq = ~ cos2x + g log sin x, g > 0, < x < tx. (2.13) 

Again by the similarity transformation in terms of the pseudo ground state wavefunction 
A<j) = 0, 0o {x) = e w{x \ we obtain 

ii = <i>? oHofo = -JfL _ 2 ^MA_ + AaM sin2 (214) 

ax" 1 ax ax 

+ (2a sin 2x - 2g cot x) — + 4aA4 sin 2 x. (2.15) 



doc doc 
Needless to say 0o has no node or singularity and it is square integrable J* (${x)dx < oo. 
The action of the Hamiltonian 7i (12.141) on monomials of sin a; reads 

H sin n x (2.16) 
-n(n — 1 + 2p) sin n ~ 2 x + n(n + 2g + 4a) sin n x + Aa(M - n) sin n+2 x, n < M - 2, 
-M{M-1 + 2g) sm M ~ 2 x + M(M + 2g + 4a) sin^x, n = M. 

Since the parity is conserved, it is now obvious that 7i keeps the polynomial space Vm 
invariant, 

UV M C V M , (2.17) 

def / Span [l,sin 2 x, ...,sin 2fc x, ..^sin^x] , Al : even, 

^ _ ) o r • -3 ■ 2fc+i -mi ka l j (2.1s) 

I Span [sin x, sin x, . . . , sin T x, . . . , sin xj , Al : odd, 

/ A/1/2 + 1, Al :even, 
dimV^ = | (> j + m M ; Qdd; (2.19) 

and that 7i is again a tri-diagonal matrix. Thus we can obtain a finite number (dimVx) of 
exact eigenvalues and eigenfunctions in the same way as in the sextic potential Hamiltonian 
f[2~Tj) case. 

2.2 'Discrete' Quantum Mechanics 

2.2.1 Difference Equation Analogue of Harmonic Oscillator Deformed by Sextic 
Potential I 

A difference analogue of the sextic potential Hamiltonian (12. ip is 

H = y / V{x) e~ id * ^fV{x)* + y/V(x)* e id * ^/VXx) - (V(x) + V(x)*) +a M {x) , (2.20) 
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= A^A + a M {x), a M {x) = 2Mx 2 , (2.21) 

^ = y/V{x) e~i 9 * - v^R*e^, A = e-^v^F-e^V^). ( 2 - 22 ) 
y(x) d = (a + ix)(b + ix)V {x), V (x) = c + ix, a,b,ceR+. (2.23) 

Here as usual V(x)* is the complex conjugate of V(x). If V is replaced by V and the last 
term in (I2.20p . a^(x), is removed, 7i becomes the exactly solvable Hamiltonian of a differ- 
ence analogue of the harmonic oscillator, or the deformed harmonic oscillator in 'discrete' 
quantum mechanics [12]. Its eigenf unctions consist of the Meixner-Pollaczek polynomial, 
which is a deformation of the Hermite polynomial [T2], [18] . The quadratic polynomial factor 
(a + ix) (b + ix) can be considered as multiplicative deformation, although the parameters a, 
b and c are on the equal footing. On the other hand one can consider it as a multiplicative 
deformation by a linear polynomial in x: 

V(x) = (a + ix)V i(x), Voi(x) = f (6 + ix)(c + ix), 

with Voi describing another difference version of an exactly solvable analogue of the harmonic 
oscillator [12]. Its eigenf unctions consist of the continuous Hahn polynomial. 

Next let us introduce the similarity transformation in terms of the pseudo ground state 
wavefunction (f>o(x): 

4>o(x) = f \/T(a + ix)T(a — ix)T{b + ix)T{b — ix)T(c + ix)T(c — ix), (2.24) 
U A ^ (j) 1 o U o O = V(x) (e~ ia ' - 1) + V(x)* {e id * - l) + 2Mx 2 . (2.25) 

It is obvious that 4>o has no node and that it is square integrable. As in the ordinary quantum 
mechanics, 4>o(x) is annihilated by the A operator (I2.22p 

= AMx) = (e-i 9 *y/V(x)* - ei 9x y/V(x)^ <j) {x). (2.26) 
It is rather trivial to verify the action of the Hamiltonian Ti. (12.251) on monomials of x: 

^ J a n>j x^-V, n < M - 2, a n>j £ R, 

Here [m] is the standard Gauss' symbol denoting the greatest integer not exceeding or equal 
to m. Since the parity is conserved, 7i keeps the polynomial space Vm invariant, 

HV M C V m , (2.28) 



6 



v def | Span [1, x 2 , . . . , x 2k , . . . , x M ] , .M : even, 
1 Span [x, x 3 , . . . , x 2fc+1 , . . . ,x M ] , M : odd, 

/ jM/2 + 1, M:even, . . 

dimV^ = | {M + 1)/2> ^ ; Qdd; (2.30) 

but 7i is no longer a tri-diagonal matrix, (H)jk ^ 0, j > k — 1. The Hamiltonian 7i ( 12.211) 
is obviously hermitian (self-adjoint) and all the eigenvalues are real and eigenfunctions can 
be chosen real. We can obtain a finite number of exact eigenvalues and eigenfunctions by 
sweeping in a similar way as in the sextic potential case (12.101) . (12. lip . The oscillation theo- 
rem linking the number of eigenvalues (from the ground state) to the zeros of eigenfunctions 
does not hold in the difference equations. The square integrability of all the eigenfunctions 
^{^dx < oo is obvious. 

2.2.2 Difference Equation Analogue of Harmonic Oscillator Deformed by Sextic 
Potential II 

Another difference analogue of the sextic potential Hamiltonian ( 12. lft has the same form as 
(12.201) . ( 12.211) and (12.22}) . with only the potential function V(x) and the compensation term 
atM( x ) are different: 

V(x) = f (a + ix)(b + ix)Vo(x), Vq(x) = f (c + ix){d + ix), a, b, c, d G R+, (2.31) 
a M (x) = M (M - 1 + 2(a + b + c + d))x 2 . (2.32) 

This Hamiltonian can be considered as a deformation by a quadratic polynomial factor (a + 
ix)(b + ix) of the exactly solvable 'discrete' quantum mechanics having the continuous Hahn 
polynomials as eigenfunctions [12j, another difference analogue of the harmonic oscillator. 
See the comments in section 5 of [21J. 

The pseudo ground state wavefunction <po(x) annihilated by the A operator A(J)q = 
reads 

0o(x) (2.33) 
= f \AX<z + ix)T(a — ix)Y{b + ix)Y{b — ix)T(c + ix)T(c — ix)T(d + ix)T(d — ix). 

Again it has no node and it is square integrable. The similarity transformed Hamiltonian 
acting on the polynomial space is 

U = 0" 1 o H o 0o = V(x) (e-^ - 1) + V(x)* (e id * - l) 

+M (M - 1 + 2(a + b + c + d))x 2 . (2.34) 



It is straightforward to verify the relationship (12.271) and to establish the existence of the 
invariant polynomial subspaces of given parity (12.281) . (12.291) and (12.301) . The hermiticity of 
the Hamiltonian and the square integrability of the eigenf unctions also hold. Thus another 
example of quasi exactly solvable difference equation is established. 

3 Other Examples 

The other two examples are the difference equation analogues of the harmonic oscillator 
with the centrifugal potential deformed by the sextic potential. There are two types corre- 
sponding to the linear and quadratic polynomial deformations. The corresponding exactly 
solvable difference equation has the Wilson polynomial [T2j [9], [10] as the eigenf unctions. The 
last example is the difference analogue of the model discussed in subsection 12.1.2] 1 / sin 2 x 
potential deformed by a cos 2x potential. In this case the corresponding exactly solvable dif- 
ference equation has the Askey- Wilson polynomials [T2l [91 [10] as eigenf unctions. The basic 
idea for showing quasi exact solvability is almost the same as shown above. 

3.1 Difference Equation Analogues of Harmonic Oscillator With 
Centrifugal Potential Deformed by Sextic Potential 

The Hamiltonians have the same form as (12.201) . (12.211) and (I2.22p . with only the potential 
function V(x) and the compensation term a^(x) are different: 

Type I: V(x) = (b + ix)V (x) 7 a M (x)=Mx 2 J (3.1) 
Type II: V{x) = (a + ix)(b + ix)V (x), (3.2) 



a M (x) = M(M-l + (a + b + c + d + e + f))x 2 , (3.3) 



with a common Vq(x) 



rr f s dcf (c + ix)(d + ix)(e + ix)(f + ix) 

V (x) = 2ix(2ix + 1) ' d,e, f £ M + — {1/2}. (3.4) 

None of the parameters a, b, c, d, e or / should take the value 1/2, since it would cancel the 
denominator. Because of the centrifugal barrier, the dynamics is constrained to a half line; 
< x < oo. The type I case can also be considered as a quadratic polynomial deformation 
of the exactly solvable dynamics with Vqi(x): 

rp r f /f \ dcf (u , • \( , • u/ ( \ t/ / \ dcf {d + ix){e + ix){f + ix) 
Type I : V(x) = {b + tx)(c + tx)V i{x), V 01 (x) = 2ix(2ix + 1) ' ( 3 - 5 ) 



which has the continuous dual Hahn polynomials [121 El CEO] as eigenfunctions. This re- 
interpretation does not change the dynamics, since the Hamiltonian and A and operators 
depend on V(x). 

The pseudo ground state wavefunction <po(x) is determined as the zero mode of the A 
operator A(j) = 0: 



Type I: <fo(x) = V , 3.6 

y/T(2ix)T(-2ix) 



Type II: fa(x) = V , 3.7 

V ; y/T(2ix)T(-2ix) 

in which the numbering of the parameters 

def def , def def , def def „ /o o\ 

ai = a, d2 = b, a% = c, CI4 = a, a§ = e, a^ = /, (,3-oj 

is used. It is obvious that both O have no node in the half line < x < oo. 

The similarity transformed Hamiltonian acting on the polynomial space has the same 
form as before (12.341) 

H d £ ( j)- 1 oHo(l) = V(x)(e- id * -l)+V(x)* (e id * -l)+a M (x)- (3.9) 

Although the potential V(x) has the harmful looking denominator l/{2ix(2ix + 1)}, it is 
straightforward to verify that 7i maps a polynomial in x 2 into another: 

~ 9 f ESvi^ 2 " 23 ', n<M-l, a nJ eR, 
<Hx 2n = { J ~° (3.10) 

This is because Vq, which has the above denominator, keeps the polynomial subspace of 
any even degree invariant, reflecting the exact solvability. This establishes that 7i keeps the 
polynomial space Vm invariant, 

HV M Q V m , (3.11) 
V M = Span [l,x 2 , . . . ,x 2k , . . . ,x 2M ] , dimV^ = M + 1. (3.12) 

The hermiticity of the Hamiltonians is obvious and the square-integrability of the eigenfunc- 
tions J °° (f) 2 (x)dx < oo holds true. This establishes the quasi exact solvability. 
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The corresponding quantum mechanical system has the prepotential W(x) and the com- 
pensation term ajn(s) as 

W(x) = — ^x 4 — ^x 2 + glogx, a, 6, g > 0, < x < oo, (3.13) 
om^) = -glMx 2 . (3.14) 

This and the above two difference analogue systems share the same invariant polynomial 
subspace ( 13. lip . (I3.12p . The undeformed exactly solvable system, i.e. (13.131) with a = 0, has 
the Laguerre polynomials as eigenfunctions. The corresponding undeformed exactly solvable 
difference equations determined by V (13.41) . Vbi (I3.5P have the Wilson and the continuous 
dual Hahn polynomials as eigenfunctions. These are three and two parameter deformation 
of the Laguerre polynomial [12], [9] [10] • 

3.2 Difference Equation Analogue of 1/ sin 2 a; Potential Deformed 
by cos2x Potential 

This system is a quasi exactly solvable deformation of the exactly solvable dynamics which 
has the Askey- Wilson polynomials [T2l [9j [10] as eigenfunctions. Let us first introduce the 
variables and notation appropriate for the Askey- Wilson polynomials. We use variables 6, x 
and z, which are related as 

0<9<tt, x = cosfl, z = e ie . (3.15) 

The dynamical variable is 6 and the inner product is (f\g) = ft d9f{6)*g(6). We denote 
D d = z£. Then q D is a g-shift operator, q D f(z) = f(qz), with < q < 1. The Hamiltonian 
is obtained by deforming the potential function Vq(z) by a linear polynomial in z: 

H = ^V(z)q D ^V(zy + VVWQ~ D VW)-(V(z)+V(z)*) + a M (z), (3.16) 
= A^A + a M (z), a M {z) = -abcdeq- 1 {l-q M ){z + - z ), (3.17) 

^ d ^ f -i(yV{z)q% -v^?^), A^t(q%^vW-q~ f Vnz))^) 

xr/ \ def ( u// v T //\def (1 -bz)(l -cz)(l -dz)(l -ez) 

V(z) = {l-az)V {z), Vo(z) = — , 3.19 

(1 — z z )(l — qz l ) 

-1 < a,b,c,d,e < 1. (3.20) 
The pseudo ground state wavefunction <fio(z) is determined as the zero mode of the A 
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operator A(f) = 0: 



'■'"< = ^7 Hi h r^~m i r-i-. ( 3 - 21 ) 

(az, az L ,bz,oz L ,cz,cz L ,dz,dz 1 ,ez,ez 1 ;q) 00 

where (a 1; ■ • • , a m ; g)^ = fl^Li II^lo^ — djq n ). Obviously (fio has no node or singularity in 
< 6 < 7i. We look for exact eigenvalues and eigenfunctions of the Hamiltonian (13.161) in 
the form: 

= S(/>, <j>{z) = <^{z)P M {x), (3.22) 

in which Pj^(x) is a degree M. polynomial in x or in z + 1/ z = 2COS6 1 = 2x. The similarity 
transformed Hamiltonian acting on the polynomial space has the form 

H = O 1 o U o 0o = V{z) {q D - 1) + V{z)* {q- D - 1) - abcdeq-\l - q M \z + i). (3.23) 

Without the deformation factor 1 — and the compensation term, the above Hamiltonian 
Ti keeps the polynomial subspace in z + 1/z of any degree invariant. It is straightforward to 
show 

HV M C V m , (3.24) 



i / i \ fc / i \ ^ 



= Span 



, dimVx = .M + 1.(3.25) 



The hermiticity of the Hamiltonian is obvious and the square-integrability of the eigenfunc- 
tions holds also true. This establishes the quasi exact solvability. 

4 Summary and Comments 

First let us summarise: Five explicit examples of quasi exactly solvable difference equations 
of one degree of freedom are derived by multiplicatively deforming the the known exactly 
solvable difference equations for the Meixner-Pollaczek, continuous Hahn, continuous dual 
Hahn, Wilson and Askey- Wilson polynomials. The finite dimensional Hamiltonian matrix, 
no longer tri-diagonal, can be solved exactly by sweeping. All the eigenvalues and eigen- 
functions are real, but the oscillation theorem, connecting the excitation level to the number 
of zeros, does not hold. Similarity and contrast with the known QES examples in ordinary 
quantum mechanics; the harmonic oscillator (with or without the centrifugal potential) de- 
formed by a sextic potential and the 1/ sin 2 a; potential deformed by a cos 2x potential, are 
demonstrated in some detail. 
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A few comments are in order. First let us stress that the mere existence of the finite 
dimensional invariant polynomial subspace 7{Vm Q Vm is n °t sufficient for the quasi ex- 
act solvability. The theory must be endowed with a pseudo ground state wavefunction O > 
which must be nodeless and square integrable. Moreover, the reverse similarity transformed 
Hamiltonian 7i = f O ° ft ° 0(7 1 must be hermitian, in order to guarantee the real spectrum. 

Let us consider an additively deformed potential with a cubic term 

V(x) = ax 3 + V (x), a M (x) = f aM{M - l)x, aeR, (4.1) 

with an exactly solvable Vo, for example, 

V Q (x) = b + ix, or (b + ix)(c + ix), b, c G K+, (4.2) 

corresponding to the Meixner-Pollaczek and the continuous Hahn polynomials. It is rather 
trivial to verify the existence of the invariant polynomial subspace: 

H A £ V{x) (e- td * - 1) + V{x)* (e id * - l) + aM{M - l)x, (4.3) 
HV M C V m , V m = Span [l,x,x 2 , . . . ,x k , . . . ,x M ] . (4.4) 

Although we have not been able to derive the explicit form of the pseudo ground state 
wavefunction <f>o(x) as a solution of A<fto = 0, it seems rather unlikely that the 0o satisfies 
the above mentioned criteria. This is because the corresponding quantum mechanical case, 
the quartic oscillator, 

b 

W(x) = ax 3 - -x 2 , a M (x)=6aMx, (4.5) 

also have the invariant polynomial subspace (I4.4p . The square integrability of (j>o(x) = e w ^ 
does not hold whichever sign a might take. 

Let us discuss the hermiticity of the Hamiltonians (12.201) and (13.161) . The hermiticity 
means (g\Ttf) = (7~tg\f) for a given inner product (g\f) for arbitrary elements / and g in 
a certain dense subspace of the appropriate Hilbert space. The obvious choice for such a 
subspace is spanned by the 'pseudo ground state' wavefunction O times polynomials. The 
types of the polynomials are: 

(a) : polynomials in x for the Hamiltonians in section 12.21 (4-6) 

/oo 
g(x)*f(x)dx, f(x) = Mx)P(x), g(x)=<j> (z)Q(x), (4-7) 
-oo 

(b) : polynomials in x 2 for the Hamiltonians in section I3.1[ (4-8) 
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(g\f) = / g(x)*f(x)dx, f(x) = Mx)P(x 2 ), 9(x) = <l> {x)Q{x 2 ), (4.9) 








The Hamiltonians (12.201) and (13. 16ft consist of three parts: 



n = n 1 + n 2 + n 3 , n 3 = a M -(v + v*) 



(4.12) 



and 



Hi 



y/V(x) e~ id * y/V(x)*, H 2 = \fV(x)* e id * y/V(x) 



for (T2~20]) 



(4.13) 
(4.14) 



for (13TT6D . 



It is obvious that 7i 3 is hermitian by itself. When 7i\ acts on /, the argument is shifted from 
x to x — % or from 9 to 9 — i log q. With the compensating change of integration variable 
from x to x + % or from 9 to 9 + zlogg one can formally show (g\H\f) = (Ji-\g\f) in a 
straightforward way. Similarly we have (g\Ti.2f) = (H2g\f) by another change of integration 
variable. This is the 'formal hermiticity.' 

Actually, the shift of integration variable, to be realised by the Cauchy integral, would 
involve additional integration contours: 



It should be noted that all the singularities arising from V and V* in cases (b) and (c) 
are cancelled by the zeros coming from the pseudo ground state wavefunctions 0o an d 4>q, 
and the Cauchy integration formula applies in all cases. The contribution of the additional 
contour integrals (I4.15I) - (I4.17I) cancel with each other and the shifts of integration variables 
is justified. 

To be more specific, the contribution from the contours at infinity in (a) and (b) vanish 
identically due to the strong damping by </>o and <J)q. The contribution from the two vertical 
contours in (6) passing the origin cancel with each other due to the evenness of 0o, (po(x) = 
0o(— x) and the polynomials in x 2 , P((—x) 2 ) = P(x 2 ). Likewise the contribution from the 
two vertical contours in (c) passing the origin cancel with each other. Those in (c) passing 



(a) : (— oo, ±i — oo), (+oo, ±i + oo) for the Hamiltonians in section I2~2l (4.15) 
(6) : (0, ±i), (+oo, ±z + oo) for the Hamiltonians in section [3TT1 (4.16) 

(c) : (0, ±i logg), (ir, n ± i log q) for the Hamiltonians in section I3T2"1 (4-17) 
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7r also cancel with each other due to the 2-7T periodicity of <po, 4>o{@) — 0o(2vr + 9). Note that 
the hermiticity in (b) and (c) cases holds only for the sum Tii + 7i 2 . This concludes the 
comments on hermiticity. 

It is now evident that the scope of the present method, deformation of W or V for gener- 
ating a quasi exactly solvable system from an exactly solvable dynamics, is rather limited. It 
is highly unlikely to get a quasi exactly solvable system, if W contains a term higher than x 5 
or cos4x, or if V has a form rij=i( a j + n > 5 or (jYj=i( a j + / {2ix(2ix + l)}, n > 7. 
However, generalisation to multi-particle difference equations, the Ruijsenaars-Schneider-van 
Diejen systems [HE], is possible [2TJ. The corresponding ordinary quantum mechanical sys- 
tems, i.e. quasi-exactly solvable Calogero- Sutherland systems [21 E] are derived by Sasaki 
and Takasaki [T5] . 

For simplicity of presentation, we have restricted the parameters a, b, c, d, e and / in V(x) 
to be real, positive etc. In most cases the ranges of parameters can be relaxed without losing 
quasi exact solvability. For example, in the first example of subsection 12.2.11 a configuration 
a > and b and c are complex conjugate, b = c* with positive real parts, is also possible. 

It is a challenge to see if these newly derived quasi exactly solvable difference equations can 
be understood by the existing ideas of QES; the s/(2) algebra [H] and/or its deformations, 
various generalised super symmetry ideas [19] . Can we use the Bethe ansatz method to solve 
these problems? Are these QES systems equivalent to some spin systems [20] ? 
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